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to be in accord with the pyramid break-up pattern. 
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Exploring the correlated electronic motion during ion- 
ization of multi-electron atoms and molecules, for ener- 
gies close to the ionization threshold, is a fundamental 
and challenging task in physics. This electronic corre- 
lation has been a topic of intense interest, for most re- 
cent see ref. [H since Wannier's pioneering work Q. 
According to Wannier's law a oc YP for excess energies 
E — > 0, where a is the cross-section of the process in- 
volved and /? depends on the number and type of par- 
ticles involved in the break-up process. Using classical 
mechanics, Wannier also showed that two electrons mov- 
ing in the Coulomb field of an ion escape back-to-back for 
energies E^O. Extending Wannier's work, later studies 
predicted a three-electron break-up on the vertices of an 
equilateral triangle [H, H[ and a four-electron break-up on 
the vertices of a tetrahedron |1, 0] . While these highest- 
symmetry break-up patterns were predicted for E — > 
it is generally expected that they also prevail for excess 
energies a few eV above threshold where the threshold 
Wannier exponent [3 is still retrieved. In this work we 
show that this is not true. 

We show that for single-photon quadruple ionization 
(QI) from the ground state of Be the prevailing break-up 
pattern a few eV above threshold is different than the one 
predicted for E — > 0. For single-photon triple ionization 
from the ground state of Li we have already shown that 
the break-up pattern is not the expected "triangle" but a 
T-shape a few eV above threshold [7[ . In the T-shape two 
electrons escape back-to-back while the third electron es- 
capes at 90° with respect to the other two electrons. Very 
recently, further evidence for the T-shape was provided 
by fully quantum mechanical calculations for energies 5 
eV above the triple ionization threshold of the ground 
state of Li Q . The previously predicted "triangle" pat- 
tern was, however, observed in recent (e, 3e) coincidence 
measurements for electron-impact on the ground state of 
He Q . The above reinforce a prediction we made in (Io| 
that the three-electron break-up pattern depends on the 
initial state and can be either a T-shape or a "triangle" . 

In the current work, we present evidence that for 



single-photon QI from the ground state of Be, a few 
eV above threshold, the prevailing break-up pattern is 
a pyramid. That is, the three electrons escape on the 
vertices of an equilateral triangle at 120° from each other 
and the other electron escapes perpendicular to the plane 
of the three electrons. Our prediction differs from the 
symmetric four-electron escape on the vertices of a tetra- 
hedron predicted in the limit E — > @, @]. Moreover, 
uncovering the physical mechanisms of QI, we express 
the multi-electron escape dynamics in terms of momen- 
tum transferring attosecond collision sequences. Thus, 
besides addressing a fundamental law of physics, we also 
elucidate correlated electronic motion in multi-electron 
escape. This is of high interest since the advent of ultra- 
short and intense laser pulses has brought time-resolving 
correlated electron dynamics in intra-atomic ionization 
processes at the forefront of Attosecond Science llHl5|. 



Given computational capabilities, addressing four- 
electron escape with quantum mechanical techniques is 
currently out of reach [ljjj . Classical mechanics is jus- 
tified for excess energies close to threshold as detailed 
in the original work of Wannier 0] and in subsequent 
work on two 17 1 and three-electron escape 0, Qij . We 
tackle quadruple photoionization using the quasiclassi- 
cal technique — quasi due to the choice of initial state — 
detailed in [l8| . Briefly using the Classical Trajectory 
Monte Carlo method HI, 2(j, we propagate in time the 
full five-body Hamiltonian accounting for all interactions 
to all orders. In addition, we use a Wigner 2l| distribu- 
tion for setting-up the initial phase space of the bound 
electrons. We compute the probability for QI, P 4+ , for 
excess energies ranging from 3 cV to 10 eV. 3 eV is close 
to threshold and the computational time involved is not 
prohibitive for obtaining good statistics. 10 cV is an up- 
per bound estimate of excess energies where the Wannier 
exponent (3 can still be retrieved. Using our data for P 4+ 
from 3 eV to 10 eV in steps of 1 eV we find (3 equal to 
94% of the theoretically predicted value of 3.288 |5|. In 
the framework of Wannier's theory, in what follows we 
discuss our results for 3 eV and 10 eV. 



2 




(degrees) 

FIG. 1. Probability for two electrons to escape with an inter- 
electronic angle 9 for excess energies 3 eV (black dots) and 10 
eV (blue dots). To guide the eye, for each excess energy, we 
connect the dots representing our data with a line. 

To identify the four-electron escape pattern we focus 
on an observable that naturally encompasses electronic 
correlation. Such an observable is the probability for two 
electrons to escape with an inter-electronic angle 9 — we 
refer to it as angular correlation probability C(6). In 
Fig. [TJ we plot C(B) for 3 eV and 10 eV excess energies. 
Given that P 4+ is 1.8xlCT 10 for 3 eV and 7.3xl(T 9 for 
10 eV the computational task involved is immense. Nev- 
ertheless, to provide good accuracy, for each excess en- 
ergy we consider, our results involve roughly 10 4 quadru- 
ple ionization events. In Fig.[TJ we see that for 10 eV C(9) 
has two peaks: one around 74° and a second one around 
100°-125°. However, for 3 eV it is not clear whether only 
one or two less — compared to 10 eV — pronounced peaks 
are present in the range 80°-112°. In Fig.[TJ C(0) is plot- 
ted using 28 bins for 8. We choose the bin size so that 
the double peak structure in C(9) is best resolved given 
the limitations imposed by our statistics. 

To what four-electron escape geometry does the shape 
of C(9) correspond to? A tetrahedron pattern with all 
electrons escaping at 109.5° from each other would result 
in a single peak in C(9). A pyramid pattern with three 
electrons escaping at 120° from each other and the other 
electron escaping at 90° from the three electrons would 
result in two peaks in C(9). Hence, the double peak in 
G(9) (Fig.[T]) for 10 eV is more consistent with a pyramid- 
shape while for 3 eV the shape of C(0) does not provide 
conclusive evidence for the prevailing escape geometry. 

We next elucidate the physical mechanisms of QI and 
provide conclusive evidence for the break-up patterns the 
four electrons follow. How does the photo-electron re- 
distribute the energy it gains from the photon to the 
remaining three electrons? This is a natural question 
in the framework of classical mechanics where the elec- 
trons undergo soft collisions mediated by Coulomb forces. 
Does redistribution of energy take place through one si- 
multaneous collision between all 4 electrons or through a 
sequence of collisions? To answer this question, we use 



a classification scheme similar to the one we first intro- 
duced in the context of three-electron escape following 
single-photon absorption from the ground state of Li Q . 
That is, we define a collision between electrons i and j — 
labeling it as ij — through the momentum transfer 

Du = r 2 w( rij )dt (i) 

Jtl 

under the condition that V(rjj(tk)) are local min- 
ima in time with t2 > ti while ry = In — rj | and 
V(rjj) = l/|ri — rj|. During the time interval ti < t < t2 
all five particles interact with each other. Hence, the 
above definition is meaningful if the collision redistributes 
energy primarily within the three-body subsystem that 
includes the nucleus and the electrons i and j. For au- 
tomated identification of the collisions, we need sensitiv- 
ity thresholds to register only the important collisions 
for the quadruple events. We do so for each individual 
QI trajectory by forming the maximum D = maxi^j|Dy| 
and registering only collisions with |Dy|/D > 5 where 
i,j = 1,2,3,4. We introduce another sensitivity thresh- 
old for how "sharp" a collision is. Namely, if elec- 
tron i gains energy through more than one collisions, we 
find the maximum AV; = max i/j (V(r ij ) max - V(r i j) min ), 
with V(r ij ) max / min the max/min value of V(ry(t)) 
for ti < t < t2, and register only collisions satisfying 
(V(rij) max - V(r;j) min )/AVi > tfi. Wc have checked that 
our results and conclusions do not change for different 
values of S and 5\) we choose 8 = 1/12 and 5i = 1/8. 

According to this classification scheme we find that 
electrons 2, 3 and 4 gain sufficient energy to leave the 
atom through two prevailing ionization routes. We de- 
note by electron 1 the photo-electron (from an Is orbital), 
by 2 the other Is electron and by 3 and 4 the two 2s elec- 
trons. In the first route the photo-electron 1 knocks-out 
first electron 2 and then proceeds to knock-out electrons 

3 and 4. That is, first a collision 12 takes place very 
early in time and roughly 24 attoscconds later collisions 
13 and 14 occur. With collisions 13 and 14 taking place 
close in time we find that a fourth collision 34 can oc- 
cur in addition to the previous 3 collisions. We refer to 
this ionization route where the photo-electron transfers 
energy to both electrons 3 and 4 as Si = {12, 13, 14}. In 
the second route, the photo-electron 1 first knocks-out 
electron 2 through the collision 12. Then, electron 2 be- 
comes the new impacting electron knocking-out, roughly 
24 attoscconds later, electrons 3 and 4 through the colli- 
sions 23 and 24. With collisions 23 and 24 taking place 
close in time a fourth collision 34 can occur in addition to 
the previous 3 collisions. We refer to this ionization route 
where electron 2 transfers energy to both electrons 3 and 

4 as S2 = {12,23,24}. si accounts for 41% and S2 for 
24/26% of all QI events for 3/10 eV. Using this scheme 
of momentum transferring attosecond collision sequences 
we have thus obtained a physical picture of the correlated 
electronic motion in an intra-atomic ionization process. 
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FIG. 2. Same as Fig. [1] but for each inter-electronic pair 8ij for the ionization routes Si (top row) and S2 (bottom row). 



Further, this scheme offers insight in choosing the ap- 
propriate asymptotic observables for inferring the tem- 
poral profile of electron-electron collision dynamics [22| . 
This is important since developing pump-probe schemes 
to time-resolve correlated multi-electron escape is one of 
the current challenges facing Attoscience j23|. 

In what direction do the four electrons escape in path- 
ways Si and S2? For Si, at the time when all electrons to 
be ionized have received enough energy to leave the atom 
the spatial electron distribution, we refer to it as transient 
threshold configuration TTC (lOj, is ri pa r3 pa r 4 7^ r2. 
That is, the last colliding electrons 1, 3 and 4 have 
ri pa i'3 pa r 4 which is close to the fixed point (see below) 
of the four-body Coulomb problem — three electrons and 
the nucleus. Thus, one expects that electrons 1, 3 and 4 
will escape symmetrically on a plane at 120° from each 
other. In Fig. [5] (top row) we plot C(9) for each of the six 
inter-electronic angles of escape using only the QI events 
corresponding to the Si, i.e, we plot C Bl (9). Indeed, we 
see that C Sl (0) for #13, 6*14 and #34 peak around 115°, 
both for 3 eV and 10 eV, corresponding to electrons 1, 
3 and 4 escaping on the vertices of a "triangle". (We 
note that the distributions in Fig. [T] and Fig. [5] are con- 
voluted by the polar angle volume element sin 9 resulting 
in a peak at 120° being shifted to slightly smaller an- 
gles while a peak at 90° is not affected). In addition, wc 
see that C S1 (0) for 12 , 9 2 3 and 6* 2 4 peak around 65°-75° 
and 75°-85° for 10 cV and 3 eV, respectively. Note that 
the shifting of the peak at smaller angles from 65°-75° 
for 10 eV to 75°-85° for 3 eV shows a tendency towards 
the pyramid-consistent angle of 90°. Thus, the distribu- 
tions in Fig. [5] (top row) for the Si ionization route are 
consistent with the pyramid-shape shown in Fig. [3] (a) . 
Similarly for the ionization route S2, C S2 (9) for #23 ; $24 
and #34 peak around 115° while C S2 (9) for 612, #13 an d 
6*i4 peak around 65°-75° and 85° for 10 eV and 3 eV, re- 
spectively (Fig. [2] bottom row). These distributions are 
consistent with the pyramid-shape shown in Fig. [3] (b). 
Therefore, for the majority (65%) of QI events the four 



electrons escape on the vertices of a pyramid. 





FIG. 3. The pyramid escape geometry for four electrons cor- 
responding to collision sequences Si (a) and S2 (b). 

We now provide further evidence that if the three elec- 
trons escape on a plane at 120° from each other then the 
preferred escape geometry of the other electron is perpen- 
dicular to this plane. We do so analytically by expressing 
the five-body Hamiltonian in hypcrsphcrical coordinates 



p 2 A 2 

2 2R 2 R 



c(n) 



(2) 



where J! = (ax, a 2 , a 3 , Ox, 9-2, 9 3 , 4 , xi, X2, X3, X*) con- 
tains all angular variables describing the positions of the 
electrons and A is a function of fi and all conjugate mo- 
menta. The total Coulomb interaction V = C /R acquires 
in this form simply an angular dependent charge C(f2). 
The hypcrsphcrical coordinates are given by 



ax = arctan(^) 
a 2 = arctan( 



Xi 

X2 



q 3 = arctan(^^li) 



X3 = <?2 -1 

Xa = <Pi + < 



(3) 

where </h, 9, are the azimuthal and polar angles of the ith 
electron. Focusing on si, the TTC is n pa 1-3 pa 1-4 7^ ra. 
For simplicity we assume ri pa r3 pa r 4 << r2 resulting in 
a 3 pa (the opposite case would lead to the same result). 
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We then expand C(D) in powers of 0:3 

3 



(4) 



5. Thus, the problem of finding a stable configuration 
is that of the three-electron problem with the solution 
a\ = tt/4, a* 2 = arctan(v / 2), X* = 27r/3, \* 2 = 4tt/3 and 
6\ = #3 = 64 = 90° [10| . These values minimize C2 for 
any value # 2 . Minimizing C3 with respect to 9 2 we find 
the stable solution 9 2 = 0° which indeed corresponds to a 
The lowest order term in a 3 , is the potential term pyramid break-up geometry, which is of lower symmetry 
of the four-body Coulomb problem with Z = 4, see Eq. than a tetrahedron. 

I 



Z 



1 



sin(a 2 )sin(Q;i) sin(a 2 ) cos(ai) cos(a 2 ) sin(a 2 ) y/l - sin(2ai) (sin 9 1 sin 9 3 cos(xi) + cos 6»i cos 9 3 ) 

1 

\J sin 2 (a 2 ) sin 2 (ai) + cos 2 (a 2 ) — sin(2a 2 ) sin(ai) (sin6*i sin #4 cos(x2) + cos£q cos 6 4) 



+ 



sin 2 (a 2 ) cos 2 (ai) + cos 2 (a 2 ) — sin(2a 2 ) cos(ai) (sin #3 sin #4 cos(x 2 — xi) + cos #3 cos #4 



r 



(5) 
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FIG. 4. C(9) for the ionization routes S3 (a) and si + S2 (b). 
The lower statistics in (a) compared to (b) dictate using 18 
(a) instead of 28 (b) bins. 

It is now clear that for the Si + s 2 labeled QI events 
the four electrons escape on the vertices of a pyramid, 
see Fig. [2] Since these events account for roughly 65% of 
all QI events the pyramid-sha,pe prevails for 3 eV and 10 
eV excess energy. Why then is the double peak in C(9) 
more pronounced for 10 eV than for 3 eV in Fig.[T]? One 
reason is the following: C Sl (9) for 6*i 2 , 9 23 and # 2 4 (Fig. [2 
top row), and C S2 (9) for 6*i 2 , 9 13 and #14 (Fig. ® bottom 
row) peak at 85° for 3 eV while only at 65-75° for 10 
eV. In addition, C Sl (0) for 9 13 , 9 U and 6» 34 (Fig. [U top 
row), and C S2 (9) for 23 , # 2 4 and #34 (Fig. [5J bottom row) 
peak roughly at 115° for 3 eV and for 10 eV. Thus, by 
summing up to obtain C Sl+S2 (9), see Fig. |U(b), the two 
peaks are closer for 3 eV (at 85° and 115°) than for 10 eV 
(at 65°-75° and 115°) resulting in a stronger overlap and 
a less pronounced double peak for 3 eV. The same effect 
is also present when all QI events are considered in C(9) 
in Fig. [1] Another reason is an ionization route which 
involves at least four distinct collisions: one collision is 
12 while two of them involve electron 3 and/or 4 each 
gaining energy by both electrons 1 and 2 — we label this 
route as S3. For S3 TTC is ri ps t 3 w 14 w r 2 . This spa- 



tial distribution is close to the fixed point of the five-body 
Coulomb problem corresponding to all four electrons es- 
caping on the vertices of a tetrahedron at 109.5° from 
each other. Indeed, in Fig. @] (a) we find that C S3 (9) has 
a single peak consistent with a tetrahedron geometry. As 
expected this single peak becomes sharper with decreas- 
ing excess energy; compare C S3 (9) for 10 eV and 3eV in 
Fig. |4](a). Thus, when all ionization routes are consid- 
ered the contribution of C S3 (9) for 10 eV docs not smear 
out the double peak of C Sl+S2 (9), see Fig. |U(b), while it 
does so for 3 eV. Further contributing to the difference 
in the shape of C(9) between 3 eV and 10 eV is that the 
% contribution of S3 to all QI events increases with de- 
creasing excess energy from 7% for 10 eV to 11% for 3 
eV. It is interesting to note that while the tetrahedron 
does not prevail in the break-up geometry, as generally 
expected, it is nevertheless present. 

In conclusion, we have shown that a pyramid is the 
prevailing break-up pattern for QI by single-photon ab- 
sorption from the ground state of Be for excess energies as 
low as 3 eV above threshold. This pattern can be verified 
by future quantum mechanical and experimental studies 
of differential cross sections. Such studies have already 
been performed for three-electron atoms, see for example 
B l9l Upl . From our previous results on triple ioniza- 
tion [10( and our current on QI we conjecture that the 
four-electron break-up pattern is also initial state depen- 
dent. That is, a tetrahedron will be the break-up pattern 
for initial states where 3 electrons occupy orbitals with 
similar spatial distribution. Interestingly, we have shown 
that the pyramid-shape is in accord with a classification 
scheme of momentum transferring attosecond collisions. 
This scheme was proven to elucidate correlated multi- 
electron escape in intra-atomic processes, see also 0, H3|- 
In doing so it allows for a new perspective on how to 
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time-resolve collisional multi-electron escape dynamics, 
a problem of intense interest in Attoscience. 
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